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Résumé (en arabe) 
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Abstract 


In this thesis, We investigate the processes of eigenvalues and eigenvec- 
tors of a symmetric matrix valued process X t , where X t is the solution of a 
general SDE driven by a G —Brownian motion matrix. Stochastic differen- 
tial équations of these processes are given. This extends results obtained 
by P. Graczyk and J. Malecki (2013) [16]. 

Keywords : G—Brownian motion matrix, G—stochastic differential 
équations, random matrices, eigenvalues, eigenvectors. 
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Résumé 


Dans cette thèse, on examine les processus de valeurs propres et des 
vecteurs propres d’un processus matriciel symétrique X t . où X t est la solu¬ 
tion d’une équation différentielle stochastique générale définie à partir d’un 
G—mouvement Brownien matrice. On donne les équations différentielles 
stochastiques de ces processus. Ces travaux constituent une généralisation 
des résultats obtenus par P. Graczyk et J. Malecki (2013) (voir [16]). 

Mots-clés : G—mouvement Brownien matriciel, G—équations différen¬ 
tielles stochastiques, matrices aléatoires, valeurs propres, vecteurs propres. 
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Introduction 


La théorie des matrices aléatoires est un domaine de recherche très actif 
en mathématique moderne avec de nombreuses applications en physique 
mathématique et théorique, en analyse mathématique et également en pro¬ 
babilité. Les applications de la théorie des matrices aléatoires en statistique 
mathématique trouvent leurs origines dans les travaux de Wishart obtenus 
en 1928 sur les matrices de corrélation (c./. [33,47]). Le véritable début 
de la recherche dans le domaine des matrices aléatoires est principalement 
dû aux travaux très influents de Eugene Wigner durant les années 1950, 
motivés par des applications en physique nucléaire (voir [1, 6, 8,12, 33]). 

Plus récemment en 2013, Graczyk et Malecki [16,17] ont donné, dans 
un contexte général, le système des équations différentielles stochastiques 
(EDS) satisfaites par les valeurs propres et les vecteurs propres de la solu¬ 
tion X t , dans l’espace des matrices symétriques d’ordre n x n, d’une EDS 
gouvernée par un mouvement Brownien matrice classique de dimension 
n x n. Sous certaines conditions sur l’équation différentielle stochastique 
satisfaite par X t , ils ont établi l’existence et l’unicité de la solution des 
équations différentielles stochastiques des valeurs propres et des vecteurs 
propres et ont montré que les valeurs propres ne se rencontrent jamais. 

D’autre part, au cours de la dernière décennie, la théorie et la méthodo¬ 
logie d’espérance non linéaire ont été bien développées et ont reçu beaucoup 
d’attention dans certains domaines d’application tels que la finance, la me¬ 
sure des risques et le contrôle. Un exemple typique d’espérance non linéaire, 
appelée G—espérance a été introduit par Peng (voir [38,39,41,42,44,46]). 
Dans ce cadre de la G—espérance, un nouveau type de mouvement Brow¬ 
nien, appelé G—mouvement Brownien a été construit et le calcul stochas¬ 
tique connexe a été établi. 

Le but de cette thèse est d’unifier la notion de matrices aléatoires et 
le G—calcul stochastique pour étudier les équations différentielles stochas¬ 
tiques des valeurs propres et des vecteurs propres pour un processus matri- 
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ciel. A cet effet, on considère la G— EDS générale suivante 

dX t = g (X t ) dB t h (X t ) + h (X t ) dBjg (X t ) + a (X t ) dt + c (X t ) d (B) t 

où {B t ) t>0 est un G—mouvement Brownien matriciel de dimension n x n et 
Bj désigne la transposée de la matrice B t . 

Le processus stochastique matriciel X t prend des valeurs dans l’espace 
des matrices symétriques d’ordre n x n et les fonctions g, h, a, c : M —■> M 
agissent sur le spectre de X t . Les principales difficultés résident dans le fait 
que la G—espérance n’est pas linéaire et que (B) n’est pas un processus 
déterministe. La notion d’indépendance des variables aléatoires par rapport 
à une espérance non linéaire étant délicate, on supposera que ( B 1 B B kl ) = 
SikSjiD pour certains processus croissants bB Comme dans [16], on dérive un 
système d’EDS pour les valeurs propres et les vecteurs propres de la solution 
de X t , qui est garantie par les conditions de Lipschitz et de croissance 
linéaires, et on prouve que les valeurs propres ne se rencontrent jamais. 

Le reste de la thèse est rédigé sous forme de trois chapitres, d’un ap¬ 
pendice et d’une conclusion générale et est organisé comme suit. 

Dans le chapitre 1, on rappelle le cadre de la G—espérance et on adapte 
ce concept en fonction de notre objectif. En outre, on donne les propriétés 
connexes du G—mouvement Brownien matriciel et la G—formule d’Itô. 

Dans le chapitre 2, on donne nos principaux résultats à savoir les équa¬ 
tions différentielles stochastiques des valeurs propres et des vecteurs propres 
de X, ainsi que le fait que le temps de collision est infini quasi sûrement. 

Dans le chapitre 3, on énonce et on démontre le théorème d’existence 
et d’unicité des solutions pour les équations différentielles stochastiques 
gouvernées par le G—mouvement Brownien matriciel. 

Dans l’appendice, on présente quelques notations d’algèbre matricielle 
qui seront utilisées tout au long de ce travail. On donne également quelques 
notions de base sur le calcul matriciel et les processus matrices en général. 

Enfin, on termine par une conclusion et quelques perspectives. 
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Sara Stihi, to appear in "Ukrainian Mathematical Journal". 
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Chapitre 1 

G*—espérance et G*—mouvement 
Brownien 


Le but de ce chapitre est de rappeler quelques résultats sur la G—espérance 
et le G—mouvement Brownien. Pour plus des détails concernant la G—espér¬ 
ance, ses applications en finance et les connaissances des mesures des risques 
cohérents, on renvoi le lecteur aux références [2, 7,13,15,23,46, 52, 56, 57]. 


1.1 G—espérance 

Soient fl un ensemble donné et 7 H. un espace linéaire des fonctions à 
valeurs réelles définies sur fl tel que : 1 G PL et X G 7 H. pour tout X G PL. 

Définition 1.1.1 Une espérance sous linéaire E sur 7 H. est une fonction 
E : PL —>• M satisfaisant les propriétés suivantes : Pour tout X, Y G PL, on 
a 

1. Monotonicité : E [X] > E [Y] si X > Y. 

2. Préservation des constantes : E [c] = c, pour tout c G R. 

3. Sous-additivité : E [X] — E [X] < E [X — Y], 

4. Homogénéité positive : E [AA"] = A E [X], pour tout A > 0. 

le triplet (fl, 77, E) est appelé espace d’espérance sous-linéaire. 

Remarque 1.1.1 L’espace Pi peut être considéré comme l’espace des va¬ 
riables aléatoires. 
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Dans ce qui suit, on considère le type d’espace des variables aléatoires 
H avec la propriété supplémentaire de stabilité par rapport aux fonctions 
bornées Lipschitziennes. Plus précisément, on suppose que si X, G 'H. i — 
1 , ...,d alors 

T (Ad, ...,X d ) E H pour tout iç E C b . Lip (R d ), 

où C b . Lip (R d ) désigne l’espace linéaire des fonctions bornées Lipschtiziennes 
sur 

Théorème 1.1.1 (voir [7,20,36],) Soit (12,7 ï.E) un espace d’espérance 
sous-linéaire. Alors il existe une famille de mesures de probabilités V sur 
(12, B (12)) telle que 

E [X] = sup E P [X] pour X e H, 

P£V 

où B (12) désigne la tribu borélienne sur 12 et Ep désigne l’espérance linéaire 
sous P. 

Naturellement, on peut définir une capacité de Choquet sur H par 
c (A) := sup P [A], A E B (12). 

P GP 

Pour plus de détails, le lecteur pourra consulter [27, 36]. 

Définition 1.1.2 On dit qu’un ensemble A E B (12) est polaire si et seule¬ 
ment si c (A) = 0 et qu’une propriété a lieu quasi-sûrement (q.s. en abrégé) 
si elle a lieu en dehors d’un ensemble polaire. 

Remarque 1.1.2 Une propriété est dite vraie quasi-sûrement (q.s.) si et 
seulement si elle est vraie presque sûrement pour chaque P E V. 

Définition 1.1.3 Soit X = (X l5 ...,X d ) un vecteur aléatoire , où X, E H, 
i — 1,..., d. La distribution de X est donnée par la fonction F\ [.] suivante : 
pour tout E Cb.Lip (R d ), 

FxM :=E[<p (X)] 

De plus, si X' est un autre vecteur aléatoire d— dimensionnel et pour tout 
T e C b . Lip (R d ), F x [<p] = F x , [<p], X et X' sont dit identiquement distribués 

(X = X'). 

Soit Y := (Y\, oùY i E'H,i = l,...,m. Le vecteur Y est dit 

indépendant de X si pour chaque G C b j Jip (M d+m ) 

E [ v (X, y)] = E[E[ip (x, Y) U J . 

On remarque que ces nouvelles définitions sont compatibles avec celles 
du cadre classique si l'espérance E est linéaire. 
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1.1.1 Distribution G—normale et G—mouvement Brow¬ 
nien 

Après les définitions de base ci-dessus on introduit maintenant la notion 
de la distribution G—normale. 

Définition 1.1.4 Un vecteur aléatoire à d—dimension X = (Ab, X d ) 
défini sur un espace d’espérance sous-linéaire (fl, TL. E) est dit G—normalement 
distribué si pour tout a, b > 0 : 

aX + bX = Va 2 + b 2 X, 
où X est une copie indépendante de X, et 

G (A) := l -E[(AX,X)} 

ici désigne l’ensemble des matrices symétriques d’ordre d. 

Dans [40,41,45], Peng montre que A" = (Ab,..., Ab) est G—normalement 
distribué si et seulement si u ( t , x ) := E [<p (x + y/tX^\ , (t, x ) G [0, oo) xR d , 

T e C l Lip (M d ) est Punique solution de viscosité de l’équation parabolique 
aux dérivées partielles, appelée G—équation de la chaleur, suivante : 

f d t u(t,x) — G (Du(t,x )), {t, x) G [0, oo) x R d . . 

\ n (0. .r) — p (x) ’ ' ‘ ' 

où Du ( t , x) est la matrice hessienne de u {t, x ). 

G (.) : Srf —■> M est une fonction monotone et sous-linéaire sur § d , à 
partir de laquelle on peut déduire qu’il existe un sous-ensemble S G SJ 
fermé, borné et convexe tel que 

G (A) = - sup tr (AB). 

2 se e 

On écrit X ~ J\f (0; S). 

Remarque 1.1.3 Le cas réel (d = 1) correspond à S = [a 2 ;W 2 ] avec 0 < 
a<a et G = Gâ,a étant la fonction sous-linéaire paramétrée par a et et : 

G (a) = ^ (cr 2 a + — a 2 a~) , a G M, 

oùct 2 = E [A 2 ], a 2 = — E [—A" 2 ], a + = max {a, 0} et a~ = max {—a, 0} . 
Dans ce cas on écrit X ~ A/"(0; [çr 2 ;cf 2 ]) (c.f. [20], [50], [52] 
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Définition 1.1.5 Un processus (Bfi t>0 à d—dimension défini sur (12, TC, E) 
est dit G—mouvement Brownien de dimension d (réel dans le cas où d = 1) 
si les propriétés suivantes sont satisfaisantes : 

(*) B 0 = 0 ; 

(ii) pour chaque t,s > 0, l’accroissement. B t+S — B t est A/"(0; sE) —distribué 
et indépendant de (B tu ..., B tn ), pour tout n G N et 0 < t\ < ... < 
t n < t. 

Pour la simplicité de notation comme dans [3, 7, 27], on désigne par B 1 
la i eme coordonnée de B et on notera 

Bf := (a, B t ) pour tout a = (ai,..., afi T G 

Selon la définition ci-dessus, on a la proposition suivante qui est importante 
pour les développements ultérieures. 

Proposition 1.1.1 (c.f. [44]J Soit (Bfi t>0 un G—mouvement Brownien à 
d-dimension défini sur un espace d’espérance sous linéaire (12,77, E). Alors 
(Bf) t>0 est un G a —mouvement Brownien à 1—dimension, où 

Ga (a) = ^ (^Lt« + - Q_ 2 aa Ta-) 

avec 

a 2 aa r = 2G (aa T ) = E [{a, B,) 2 ] et a 2 aaT = -2G ( -aa T ) = —E [- (a, B fi 2 ] . 
En particulier, pour tout t, s > 0, Bf +S — Bf Al (0, [so_ 2 aa TiS<j 2 aaT \) . 

1.1.2 G—espérance conditionnelle 

On considère 12 = C Q (M") l’espace des fonctions continues (wfi t&R+ 
définies sur R + ,à valeurs dans M n satisfaisant lo 0 = 0, muni de la distance 
suivante : 


P pv 2 ) = £ 2- j 

i= 1 


[ max \ul — wfi ] Al 
yt£[o,i] *7 


u , ui G 12. 


Pour tout t G [0; oo), on pose 12 f := {o;. At : ce G 12}. Les espaces des fonctions 
Lipschitziennes sur 12 sont notés par : 


Lip (fi t ) = {(f (B tlM , ...,B tnM ) : ..., t n G [0; oo) , G C b . Lip (M n )} , 
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oo 

Lip (fl) = (J Lip (fl n ). 

n=l 

Ici, on utilise Gf,.^ (R”) seulement pour la commodité des techniques. En 
général, Cb.Lip (IR") peut être remplacé par l’un des espaces des fonctions 
défini sur R" suivantes. 

•L°° (R") : Espace des fonctions bornées Borel-mesurables, 

• C U nif (R") : Espace des fonctions bornées et uniformément continues, 
•Ci,Lip (R n ) : Espace des fonctions continues, localement Lipschitziennes, 
•Lip(M. n ) : Espace des fonctions Lipschitziennes sur R". 

Soit T > 0 un temps fixé. On désigne par L P G (Ht), P > 1, l’adhérence 

de Lip{Q.T ) par rapport à la norme ||X|| pG := (E [|X| p ])p . Peng dans 
[39,43,45], a construit une espérance sous-linéaire sur (fl, Lip (fl)) de telle 
sorte que le processus canonique (B t ) t>0 soit un G—mouvement Brownien 
de la manière suivante : soit (^)^ x une suite des vecteurs aléatoires à 

d —dimension sur un espace d’espérance sous-linéaire (Vl. 'H. p 'j telle que 

soit G—normalement distribuée et £ i+1 indépendante de (£ 1; ..., pour 
tout i = 1, 2 ,.... 

On introduit une espérance sous-linéaire E définie sur Lip (fl), via la 
procédure suivante : Pour tout X G Lip (fl) avec 

A = p [B tl - B to ,..., B tn - B tn _^j 

pour ip G Cb.Lip (R n ) et pour 0 = t 0 < t\ < ... < t n < oo, on pose : 


E[X\ = E (pty/ti- to£ 1, ..., y/t n - 


D’après Peng [40,43], E ainsi définie est une espérance sous-linéaire sur 
Lip (D) et (B t ) t>0 est un G—mouvement Brownien. Comme Lip( fl^) Ç 
Lip( fl), E est également une espérance sous-linéaire sur Lip{ fl T ). 


Définition 1.1.6 L’espérance conditionnelle de X par rapport, à f l tj est 
définie par 


E[X \Q tj ] ='fi{B tl ,B t2 -B tl ,...,B tj -B tj _ 1 ) (1.2) 


ou 


1p (Xi, ..., Xj) E p (xi, ..., Xj, \Jtj+ 1 tj^j- 1-1, •••’ '\/tn tn—lfil 
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Définition 1.1.7 Un vecteur aléatoire à n—dimension Y G (. Lq (fi)) ra est 
dit indépendant de O t ( t > 0), si pour tout p G Cb.Lip {Ain) on a 

E[< P {Y)\n t \ = E[<p (y)]. 

Proposition 1.1.2 Soient X, Y G L l c (fi) tels que E [Y |fi t ] = —E [-Y |fi t ] 
pour certains t G [0, T]. Alors on a 

E [X + Y |fi t ] = E [X |fi t ] + E [Y |fi t ]. 

En particulier, si E [Y |fi t ] = E [—Y |fi t ] = 0, alors 

E[x + Y\n t ] = E[x\nt]. 

Preuve, (voir [15], [20]). ■ 

Exemple 1.1.1 Pour tout a G fixé, s < t, on a 

E [B a t - B: |fi s ] = E [- (B* - B a s ) |fi s ] = 0 

E[(Bi-B‘Ÿ\a,} = a Lx(i-s), 

b[-(b«“-b;) 2 |sî s ] = 

e[(b;-b:ÿ isî.] = 3(5^7- (t — s ) 2 , 
e [- (S,“ - B»“) 4 |SÎ.] = -3d„x(i- S ) 2 . 

Exemple 1.1.2 Pour tout a G fixé, n G N, 0 < t < T, X G Lq (fi t ) et 
T e C b . Lip (K), on a 

E[Xq>{B a T -B a t )\Çl t ] = X + E[<p(BÏ-B?)\n t } + X-E[-<p(BÏ^B?)\n t ] 

= X + E [<p {B* - £“)] + X-E [-(p (B a T - S“)]. 

En particulier, on a 

E [X (B% - B?) |fi f ] = X + E [(££ - B?)} + X~E [- (. B £ - B?)} = 0. 

Ceci, combiné avec la proposition ci-dessus donne 

E [Y + X (£“ - Bî) |fi t ] = E [Y |fi t ], pour tout Y G Lq (fi). 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


On a également 


e [x ç b “ - b a t f ia; 


X + E [(£“ - B(f] 

+X~E [- (£“ - B “) 2 ] 

(r - 1 ) 


et pour tout n > 1 

E [X (B£ - B*) 211 - 1 |a] 


X + E [(££ - i?“) 2n_1 ] 
+X~E [- (B% - B() 2n ~ l ] 


\X\E (B a T _ t ) 


2n—l 


Définition 1.1.8 (voir [58, 50],) Un processus ( M t ) t>0 est appelé G—martingale 
si pour tout t G [0; T ], M t G L l G (a) G pour tout s G [0, t] on a 

E [M t |a] = M s q.s.. 


1.2 G —mouvement Brownien matrice et G— 
calcul stochastique 

Dans ce qui suit on identifie chaque matrice carrée d’ordre n à un vecteur 
de dimension n 2 et on considère maintenant 12 = Go (A4 n ) l’espace des 
fonctions continues (u t ) teR+ à valeurs dans l’ensemble des matrices carrées 
A4 n satisfaisant u 0 = 0. Les définitions du paragraphe précédent restent 
valables sur l’ensemble des matrices A4 n . 

Dans toute la suite, on considère le processus canonique (B t ) t>0 (i.e. 
B t (u>) = oj t ) défini sur (12, Lip (12), E) où 12 = C 0 (A4 n ). (B % G est donc un 
Grj o— ~mouvement Brownien réel où 

u %j i u xj 



E 


{B l \ 


ijy 


et <7; 


-E 



Soit 7r(p = {0 = tg < ti < ... < G = T} une suite de subdivisions de 
[0, T], dont le pas 

lGt) = max (K fc +i ~ t i \ '■ i = °,1 ,-,N- 1} , 
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CHAPITRE 1. G-ESPÉRANCE ET G-MOUVEMENT 

BROWNIEN 


tend vers 0 lorsque k —» oo. Par abus d’écriture on écrit i tt = {£ 0 , G, -, 
au lieu de nf = {0 = £§ < < ••• < t% = T} . Pour tout p > 1, on consi¬ 

dère l’espace des processus simples M G ° (0, T) tels que 

N -1 

M = X! M 1 [*nb+i) (*) » P° ur 0 = t 0 < ... <t N = T, 

3 =0 


V = Vt (U) = 


oùÇ j eL p G (n tj );j = 0,...,N-l. 

Pour tout p > 1, on note par M G (0, T) l’adhérence de Mff (0, T) sous 
la norme 


I e 

" T 

1 \Vt\ P dt 

J 

1 

f 

{ 

. 0 

J 


Il 't\\M q(0,T) 

Notons que M G (0, T) un espace de Banach. 

Définition 1.2.1 L’intégrale de Bochner de ( rj t ) est définie par 

T „ N-l 

/ Vt M dt = ( u ) (t j+ 1 - ij). 

0 J=° 

On donne maintenant la définition de la G—intégrale d’Itô. 

Définition 1.2.2 Pour tout a G M G ° (0, T), on définit la G—intégrale 
d’Itô par 

T r n- 1 

/ (o) : / . , 

j=0 


jV—1 

/ =2:?, (s? + , - g 

J n 


Notons que / (a) est complètement indépendante de G, et 
est une G—martingale (voir [58]). 


( t 
y 

v 0 


Oi s dBf 


0 <t 


Lemme 1.2.1 L’application I : M G ° (0, T) —*• L G (Ht) est. linéaire, conti¬ 
nue et peut-être donc prolongée par continuité à Mfi (0, T) et on a 


(13) 


E / atdBV = 0 



1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


et 



(1.4) 


Preuve. (voir [3], [25]). ■ 

Proposition 1.2.1 Soient rj,6 G Mq (0 ,T) et 0 < s < r < t < T. Alors 
on a 

(0 fVudBV = fr/udBÿ + JnJBfi, 


s s r 

t t t 

(ii) J (a? 7 u + 9 U ) dBfi = a Jr] u dB^+J9 u dB z fi , si a est bornée dans L G (fi s ), 


(ni) E 


x + frjBÿ ia 


= E [X la] pour X 6 L l G (a). 


Processus de variation quadratique 

Le processus de variation quadratique du G—mouvement Brownien est 
un processus très intéressant. On a vu que le G—mouvement Brownien est 
un processus de variance incertaine mais de moyenne certaine centrée. Cette 
incertitude est concentrée dans sa variation quadratique (IL- 1 ). Par ailleurs 
(IL 3 ) lui-même est un processus de moyenne-variance. 

Soit nfi , N — 1,2,... une suite de subdivisions de [0, t] dont le pas tend 
vers 0. On a alors 


TV—1 



Définition 1.2.3 La variation quadratique de [B\ 3 ^ t>Q est définie par la 
limite l? G (fl) suivante 


TL, 


TV—1 


Uni 


vU t) 


E ( B L 

1=0 



t 

(B?f - 2 JBÿdBÿ. 
0 
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CHAPITRE 1. G-ESPÉRANCE ET G-MOUVEMENT 

BROWNIEN 


Remarque 1.2.1 II est important de garder à l’esprit que {Bn’est pas 
un processus déterministe sauf dans le cas où afi = <Jij, alors B l] devient 
un mouvement Brownien classique. 


Proposition 1.2.2 (voir [I3])Pour tout 0 < s < t < oo, on a 

E [<B«), - (B ij ) s |SÎ,] = W? (t - s) (1.5) 

E [- «£“>, - (B 11 ),) |SÎ.] = (t - S) (1.6) 

Corollaire 1.2.1 [45] Pour tout 0 < s < t < T, on a 

4* < (< s ”> t+ . - (B 11 ),) < «i 2 t- 

Définition 1.2.4 L’intégrale d’un processus g G M G ° (0 ,T) par rapport à 
est définie par 

T „ N -1 

<?w = : = ££< - ( B ‘%) < LT > 


o 

r 1,0 


1=0 


L’application Q : M G (0 ,T) —> L l G (0 T ) est continue et peut, donc se pro¬ 
longer par continuité à Mf (0, T). 

Proposition 1.2.3 (voir [43] J Soit g G M G (0, T), alors 



' t} V" 


t 

E 

V s dB ÿ 

= E 

jm#). 


L Vo J \ 


J 

_0 


Le théorème suivant exprime les inégalités Burkholder-Davis-Gundy 
(BDG) pour les G—intégrales stochastiques relativement par rapport au 
G—mouvement brownien réel B l i (voir [13]). 

Théorème 1.2.1 1) Soient p > 1, g G Mf, (0 ,T) et 0 < s < t < T. Alors 
on a 




u 

/» 

p- 

E 

sup 

/ ^rd{B i ‘) r 



s<u<t 

J 

s 



< Gi(t-s) p 1 / E[\r) u \ p ]du 


où C i > 0 est une constante indépendante de g. 

2) Soient p > 2, g G M G (0 ,T) et 0 < s < t < T. Alors on a 

t 


( 1 . 8 ) 




U 

/» 

p- 

E 

sup 

s<u<t 

/ g r dBfi 

S 



<C2|(-s| 5 1 J E{Kf]du 

S 

où C 2 > 0 est une constante indépendante de g. 


(1.9) 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


La covariation quadratique du G— mouvement Brownien matrice 

On définit le processus de covariation quadratique par 

(B ij , B km ) := - [(B ij + B km ) - (B ij - B km >] (1.10) 


et on observe que 

N -1 


E O!, - B l) te - a 


p=0 


;y km 

b p +1 L P 


N—l 


E{[( S E + S t«) 


p=0 


B] 3 + 


pij _ r>km \ 

D tp +1 ^tp+l J 


-1 Bi j - B km 


d’où, puisque /i (ir^) -> 0, on a 


JV-1 


lim V - y) - B k ™ \ = (B i3 ,B km ) , 

V—»oo V *P+1 p / V p+1 P / ' ' 


et aussi 


p=0 


(B ij , B km ) t = - [< B ij + - (B ij - £ fcm )J 


= Bl j B km - I B\ 3 dB k 8 


B km dB l s L 


-)km _ / -okm 

s J 

0 0 

Il résulte que, pour tout 77 G Mq (0, T; M) on a 

TT T 

[n.d {B », B*”) = \j v ,d (B‘i + B k ™) - l L,d{B» - B km ) 


1.2.1 G—Formule d’Itô 

On donne maintenant la G—formule d’Itô dans un cadre général (voir 
[25,26] pour le cas vectoriel). Dans la suite, nous adoptons la convention 
de notation d’Einstein. 
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Théorème 1.2.2 Soit ip une fonction de classe C 2 sur A4 n de telle sorte 
que d x pqip, d 2 p , q , pq <p G Cb.Lip {Ain) ■ Soit X = (X 13 ) un processus d’Itô 
matriciel sur [0, T] avec 


X pq = X pq + / a pq (s) ds + / (a) d p pq (s) dB k \ (1.11) 


où a pq , 9 pq kl G Mq (0, T) et (3™ G Mf, (0, T). Alors pour chaque t G [0, T], 
on a q.s., 


t t 

r r 

-)kl 


T (X t ) - (X 0 ) = / d x „<p (X u ) Plï (u) dB k u l + / d^ip (X u ) a pq (u) du 


o 


+ / \d x pq(p (X u ) 9 pq kl ( u ) 


+\ d lv'q' x ™V ( X u) Pïj ( u ) Pki' ( u ) 


d(B ij ,B kl )l 1.12) 


Remarque 1.2.2 Notons que cette formule reste valide si X n’est, pas une 
matrice carrée. 


Le corollaire suivant est une conséquence immédiate du théorème (1.2.2). 
Dans toute la suite, on pose : 

dX pq dX™ n = P v t qij PT nkl d{ Bi PB kl ) t . 


Corollaire 1.2.2 (voir [31 ]) 

1. La formule d’intégration par partie suivante a lieu : 

d (X pq X t mn ) = dX pq X t mn + X pq dX™ n + dX^dX?", (1.13) 


2. On ad (. X pq , X™ n ) = dX pq dX™ n 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


Preuve. 1. La G— formule d’Itô appliquée avec (x, y) = xy , donne 
d(Xf q X t 


rmn\ 

J 


(i) + M AT.-v, 


V 


mn\ 




V 

av (A T, X, 


mn > 
t y 


Piï (t) PT (t) ) d (B ij , B kl ) 


dxdy /t 

= A7 n a M (t) dt + X™ n fi pq (t) dBf + X™ n e^ kl (t) d (B ij , B kl ) t 
+x t v n (t) dt + AT7r dBf + AT77 (*) <* 

= X™ n dXf q + Xï q dX™ n + fd p k q (t) p% n (t) d (B ij ,B kl ) t . 


2. Provient des égalités suivantes : 


dB\ j dB kl = d(B ij ,B kl ) t , 
dtdBfi = dtd(B ij ,B kl ) t = 0. 


Remarque 1.2.3 Dans ce qui suit, on écrit (B l] ) t au lieu de (B l fi B l Q t . 

Définition 1.2.5 Comme dans le cas classique, on définit l’intégrale de 
Stratonovich o pour deux processus d’Itô matrices X et Y : 

X o dY = XdY + \dXdY et dX o Y = dXY + ^dXdY, 

où dXdY est le produit matriciel classique. 

La proposition suivante a lieu. 

Proposition 1.2.4 (voir [51],) Soient X et Y deux processus matriciels 
d’Itô. Alors on a 

( i) La formule d’intégration par parties : 

d (XY) = XdY + dXY + dXdY (1.14) 
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1 11 


d (XY) = dX 

O Y + X O dY 

(1.15) 

dXo(YZ) = 

(dX oY)o Z 

(1.16) 

(YZ)odX = 

(Y o Z) o dX 

(1.17) 

11 

h 

ST 

o 

= dY T o X T 

(1.18) 


et 


Y T désigne la transposée de la matrice Y. 

Preuve. Soient X, Y, Z des processus matrices d’Itô. 

(i) Pour tout i , j G 1 ,n, on a 

d (XY) ij = d ( X ik Y kj J =Y^ â C Xik Y kj ) 

= ( dX ik Y kj + X ik dY kj + dX ik dY kj ) 

k 

= ( dXY) ij + (. XdY) ij + (dXdYf , 


d’où 


d (XY) = XdY + dXY + dXdY. 

(ii) Par un calcul simple, on trouve 

dX o Y + X o d,Y = dXY + -dXdY + XdY + -dXdY 

2 2 

= dXY + XdY + dXdY = d (XY). 


On a 


dXo(YZ) = dX (YZ) + -dXd (YZ) 

= dX (YZ) + \dX (dYZ + YdZ). 

D’autre part 

(dX oY)oZ = (dXY + ^dXdY^j o Z 


dX (YZ) + ^dXdYZ + ^ fdXYdZ + ^dXdYdz) 
dX (YZ) + X -dX (dYZ + YdZ ), 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 

d’où (1.16). L’égalité (1.17) se démontre de la même manière. 

La formule (1.18) provient de la définition de l’intégrale de Stratonovich 
et du fait que 

(, dXdY) T = dY T dX T . 
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Chapitre 2 

G*—ED S des valeurs propres et 
des vecteurs propres des 
processus matrices 


L’objectif de ce chapitre est de donner les équations différentielles sto¬ 
chastiques des valeurs propres et des vecteurs propres des processus qui 
sont eux même solutions d’EDS gouvernées par un G—mouvement Brow¬ 
nien matrice. Pour contourner les difficultés liées au fait que les entrées d’un 
G—mouvement Brownien matrice ne sont pas forcément indépendantes, on 
supposera que leurs covariations quadratiques sont milles. 


2.1 Le modèle 

Dans le reste de cette thèse, on suppose que B satisfait l’hypothèse 
suivante : 

(A) Il existe un processus réel et croissant V tel que (B l i,B kl Y = SikSjib J t 
q.s. pour chaque i,j, k. I G 1, n, où ô uv est le symbole de Kronecker. 

On a alors a 2 t <b{< a 2 t où a := maxcr^ et a minerjj. Notons que 

i,j i,j - 

dans le cas classique, l’hypothèse (A) est satisfaite avec lP t — t pour tout 
t G [0, T\. 

On considère la G—équation différentielle stochastique générale définie 

par 


dX t = g (X t ) dB t h (X t ) + h (X t ) dB T t g (X t ) + o (X t ) dt + c (X t ) d (B) t (2.1) 
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2.2. G-ÉQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VALEURS PROPRES 


où les fonctions g, h, a, c : M —> M opèrent sur le spectre de X t , et X 0 G S n 
l’ensemble des matrices symétriques avec n valeurs propres différentes. La 
matrice variation quadratique d (B) est définie par d (B) := dBdB. Notons 
que grâce à l’hypothèse (A), d (B) est une matrice diagonale où d (B) lJ = 
SijdbL 

Dans ce qui suit, soit § n l’ensemble des matrices symétriques de dimen¬ 
sion n identifié à M"( n + 1 )/ 2 . Soit X t la solution de l’EDS (2.1) lorsqu’elle 
existe et soit H J X t H t — A t diag(\i(t )) la forme de diagonalisation 
de X t , où IIj est une matrice orthonormale. Notons que X G S n et que 
pour tout / : 1 -> 1 on a / (I) = Hf (A) H 7 . On suppose que X 0 G § n 
et que les valeurs propres de X t soient ordonnées de manière croissante : 

Ai (t) < A 2 (t) < ... < A n (t) jusqu’au premier temps de collision défini de 
la manière suivante : 

Définition 2.1.1 On définit le premier temps de collision des valeurs propres 
par 

t = inf {t : Xi (t) = A j ( t ) pour i fi j} . 


2.2 G —équations différentielles stochastiques 
des valeurs propres 

Maintenant on est en mesure d’énoncer le résultat principal de cette 
thèse. Notons que dans notre modèle, l’équation différentielle stochastique 
à étudier (2.1) se comporte comme dans le cas linéaire [4, 5,14,16, 37]. Les 
techniques utilisées sont inspirées par le cas linéaire, où le G— mouvement 
Brownien joue le rôle d’un mouvement Brownien classique. Pour cela on 
aura besoin de la proposition suivante : 

Proposition 2.2.1 [31] Soit X t G une solution de (2.1). Alors la ma¬ 
trice des valeurs propres A t est solution de l’équation différentielle stochas¬ 
tique suivante : 

dA = dN — dA o A + A o dA, (2-2) 

où dN = H 7 o dX o H et A est la matrice logarithme stochastique de H 
définie par 

dA = H 7 o dH. (2.3) 

Preuve. Comme dans [4,16], on considère la décomposition 

A = H t XH, 
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CHAPITRE 2. G-EDS DES VALEURS PROPRES ET DES 
VECTEURS PROPRES DES PROCESSUS MATRICES 


d’où 

dA = d ( H t XH ) . (2.4) 

D’après la formule (2.4) et les notations d’Itô Stratonovich, on peut écrire 

dA = dH T o (XH) + H t o d (XH) 

= (d,H T o X) o H + H t o (dX oH + X o dH) 

= ( dH T o X) o H + H t o dX o H + H t o X o dH , 

d’où en posant dN = H T o dX o H, 

dA = ( dH T o HAH t ) oH + dN + H T o HAH T o dH. 

On obtient en utilisant les formules (1.16) et (1.17) : 

dA = dN+ dH T o HA +H T HAH T odH 
= dN+ dH T oHA + AH T odH 
= dN + dH T oHoA + AoH t odH 
= dN — dA o A + A o dA. 


Théorème 2.2.1 [51] Soit X t une solution de l’équation (2.1) telle que 
X 0 £ S n . Alors il existe un G—mouvement Brownien réel w l tel que = 
8ijV. Pour tout i,j £ 1 ,n et pour t < r le processus des valeurs propres A f 
est solution du système suivant : 

d\i = 2 g (Ai) h (Ai) Y,H ki dw k + a (Ai) dt + dV ü (2.5) 

k 

OÙ 

dV ij = 8 ijC (Ai) £ [H ki ) 2 db k + dR ij (2.6) 

k 

et 

dR ij = [<W(A fc )ù(A fc )MA,-) (2-7) 

k^j N k 

+g 2 (\ k )h(\ i )h(\ j )'EH li H lj db l 

i 

+6, j g 2 (\,)h 2 (\ t )Y.(H ,k Ÿdb l . 

I 
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STOCHASTIQUE 


On a également 


e [x ç b “ - b a t f ia; 


X + E [(£“ - B(f] 

+X~E [- (£“ - B “) 2 ] 

(r - 1 ) 


et pour tout n > 1 

E [X (B£ - B*) 211 - 1 |a] 


X + E [(££ - i?“) 2n_1 ] 
+X~E [- (B% - B() 2n ~ l ] 


\X\E (B a T _ t ) 


2n—l 


Définition 1.1.8 (voir [58, 50],) Un processus ( M t ) t>0 est appelé G—martingale 
si pour tout t G [0; T ], M t G L l G (a) G pour tout s G [0, t] on a 

E [M t |a] = M s q.s.. 


1.2 G —mouvement Brownien matrice et G— 
calcul stochastique 

Dans ce qui suit on identifie chaque matrice carrée d’ordre n à un vecteur 
de dimension n 2 et on considère maintenant 12 = Go (A4 n ) l’espace des 
fonctions continues (u t ) teR+ à valeurs dans l’ensemble des matrices carrées 
A4 n satisfaisant u 0 = 0. Les définitions du paragraphe précédent restent 
valables sur l’ensemble des matrices A4 n . 

Dans toute la suite, on considère le processus canonique (B t ) t>0 (i.e. 
B t (u>) = oj t ) défini sur (12, Lip (12), E) où 12 = C 0 (A4 n ). (B % G est donc un 
Grj o— ~mouvement Brownien réel où 

u %j i u xj 



E 


{B l \ 


ijy 


et <7; 


-E 



Soit 7r(p = {0 = tg < ti < ... < G = T} une suite de subdivisions de 
[0, T], dont le pas 

lGt) = max (K fc +i ~ t i \ '■ i = °,1 ,-,N- 1} , 
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tend vers 0 lorsque k —» oo. Par abus d’écriture on écrit i tt = {£0, G, -, 
au lieu de nf = {0 = £§ < < ••• < t% = T} . Pour tout p > 1, on consi¬ 

dère l’espace des processus simples M G ° (0, T) tels que 

N -1 

M = X! M 1 [*nb+i) (*) » P° ur 0 = t 0 < ... <t N = T, 

3 =0 


V = Vt (U) = 


oùÇ j eL p G (n tj );j = 0,...,N-l. 

Pour tout p > 1, on note par M G (0, T) l’adhérence de Mff (0, T) sous 
la norme 


I e 

" T 

1 \Vt\ P dt 

J 

1 

f 

{ 

. 0 

J 


Il 't\\M q(0,T) 

Notons que M G (0, T) un espace de Banach. 

Définition 1.2.1 L’intégrale de Bochner de ( rj t ) est définie par 

T „ N-l 

/ Vt M dt = ( u ) (t j+ 1 - ij). 

0 J=° 

On donne maintenant la définition de la G—intégrale d’Itô. 

Définition 1.2.2 Pour tout a G M G ° (0, T), on définit la G—intégrale 
d’Itô par 

T r n- 1 

/ (o) : / . , 

j=0 


jV—1 

/ =2:?, (s? + , - g 

J n 


Notons que / (a) est complètement indépendante de G, et 
est une G—martingale (voir [58]). 


( t 
y 

v 0 


Oi s dBf 


0 <t 


Lemme 1.2.1 L’application I : M G ° (0, T) —*• L G (Ht) est. linéaire, conti¬ 
nue et peut-être donc prolongée par continuité à Mfi (0, T) et on a 


( 13 ) 


E / atdBV = 0 



1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


et 



(1.4) 


Preuve. (voir [3], [25]). ■ 

Proposition 1.2.1 Soient rj,6 G Mq (0 ,T) et 0 < s < r < t < T. Alors 
on a 

(0 fVudBV = fr/udBÿ + JnJBfi, 


s s r 

t t t 

(ii) J (a? 7 u + 9 U ) dBfi = a Jr] u dB^+J9 u dB z fi , si a est bornée dans L G (fi s ), 


(ni) E 


x + frjBÿ ia 


= E [X la] pour X 6 L l G (a). 


Processus de variation quadratique 

Le processus de variation quadratique du G—mouvement Brownien est 
un processus très intéressant. On a vu que le G—mouvement Brownien est 
un processus de variance incertaine mais de moyenne certaine centrée. Cette 
incertitude est concentrée dans sa variation quadratique (IL- 1 ). Par ailleurs 
(IL 3 ) lui-même est un processus de moyenne-variance. 

Soit nfi , N — 1,2,... une suite de subdivisions de [0, t] dont le pas tend 
vers 0. On a alors 


TV—1 



Définition 1.2.3 La variation quadratique de [B\ 3 ^ t>Q est définie par la 
limite l? G (fl) suivante 


TL, 


TV—1 


Uni 


vU t) 


E ( B L 

1=0 



t 

(B?f - 2 JBÿdBÿ. 
0 
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CHAPITRE 1. G-ESPÉRANCE ET G-MOUVEMENT 

BROWNIEN 


Remarque 1.2.1 II est important de garder à l’esprit que {Bn’est pas 
un processus déterministe sauf dans le cas où afi = <Jij, alors B l] devient 
un mouvement Brownien classique. 


Proposition 1.2.2 (voir [I3])Pour tout 0 < s < t < oo, on a 

E [<B«), - (B ij ) s |SÎ,] = W? (t - s) (1.5) 

E [- «£“>, - (B 11 ),) |SÎ.] = (t - S) (1.6) 

Corollaire 1.2.1 [45] Pour tout 0 < s < t < T, on a 

4* < (< s ”> t+ . - (B 11 ),) < «i 2 t- 

Définition 1.2.4 L’intégrale d’un processus g G M G ° (0 ,T) par rapport à 
est définie par 

T „ N -1 

<?w = : = ££< - ( B ‘%) < LT > 


o 

r 1,0 


1=0 


L’application Q : M G (0 ,T) —> L l G (0 T ) est continue et peut, donc se pro¬ 
longer par continuité à Mf (0, T). 

Proposition 1.2.3 (voir [43] J Soit g G M G (0, T), alors 



' t} V" 


t 

E 

V s dB ÿ 

= E 

jm#). 


L Vo J \ 


J 

_0 


Le théorème suivant exprime les inégalités Burkholder-Davis-Gundy 
(BDG) pour les G—intégrales stochastiques relativement par rapport au 
G—mouvement brownien réel B l i (voir [13]). 

Théorème 1.2.1 1) Soient p > 1, g G Mf, (0 ,T) et 0 < s < t < T. Alors 
on a 




u 

/» 

p- 

E 

sup 

/ ^rd{B i ‘) r 



s<u<t 

J 

s 



< Gi(t-s) p 1 / E[\r) u \ p ]du 


où C i > 0 est une constante indépendante de g. 

2) Soient p > 2, g G M G (0 ,T) et 0 < s < t < T. Alors on a 

t 


( 1 . 8 ) 




U 

/» 

p- 

E 

sup 

s<u<t 

/ g r dBfi 

S 



<C2|(-s| 5 1 J E{Kf]du 

S 

où C 2 > 0 est une constante indépendante de g. 


(1.9) 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


La covariation quadratique du G— mouvement Brownien matrice 

On définit le processus de covariation quadratique par 

(B ij , B km ) := - [(B ij + B km ) - (B ij - B km >] (1.10) 


et on observe que 

N -1 


E O!, - B l) te - a 


p=0 


;y km 

b p +1 L P 


N—l 


E{[( S E + S t«) 


p=0 


B] 3 + 


pij _ r>km \ 

D tp +1 ^tp+l J 


-1 Bi j - B km 


d’où, puisque /i (ir^) -> 0, on a 


JV-1 


lim V - y) - B k ™ \ = (B i3 ,B km ) , 

V—»oo V *P+1 p / V p+1 P / ' ' 


et aussi 


p=0 


(B ij , B km ) t = - [< B ij + - (B ij - £ fcm )J 


= Bl j B km - I B\ 3 dB k 8 


B km dB l s L 


-)km _ / -okm 

s J 

0 0 

Il résulte que, pour tout 77 G Mq (0, T; M) on a 

TT T 

[n.d {B », B*”) = \j v ,d (B‘i + B k ™) - l L,d{B» - B km ) 


1.2.1 G—Formule d’Itô 

On donne maintenant la G—formule d’Itô dans un cadre général (voir 
[25,26] pour le cas vectoriel). Dans la suite, nous adoptons la convention 
de notation d’Einstein. 
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CHAPITRE 1. G-ESPÉRANCE ET G-MOUVEMENT 

BROWNIEN 


Théorème 1.2.2 Soit ip une fonction de classe C 2 sur A4 n de telle sorte 
que d x pqip, d 2 p , q , pq <p G Cb.Lip {Ain) ■ Soit X = (X 13 ) un processus d’Itô 
matriciel sur [0, T] avec 


X pq = X pq + / a pq (s) ds + / (a) d p pq (s) dB k \ (1.11) 


où a pq , 9 pq kl G Mq (0, T) et (3™ G Mf, (0, T). Alors pour chaque t G [0, T], 
on a q.s., 


t t 

r r 

-)kl 


T (X t ) - (X 0 ) = / d x „<p (X u ) Plï (u) dB k u l + / d^ip (X u ) a pq (u) du 


o 


+ / \d x pq(p (X u ) 9 pq kl ( u ) 


+\ d lv'q' x ™V ( X u) Pïj ( u ) Pki' ( u ) 


d(B ij ,B kl )l 1.12) 


Remarque 1.2.2 Notons que cette formule reste valide si X n’est, pas une 
matrice carrée. 


Le corollaire suivant est une conséquence immédiate du théorème (1.2.2). 
Dans toute la suite, on pose : 

dX pq dX™ n = P v t qij PT nkl d{ Bi PB kl ) t . 


Corollaire 1.2.2 (voir [31 ]) 

1. La formule d’intégration par partie suivante a lieu : 

d (X pq X t mn ) = dX pq X t mn + X pq dX™ n + dX^dX?", (1.13) 


2. On ad (. X pq , X™ n ) = dX pq dX™ n 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 


Preuve. 1. La G— formule d’Itô appliquée avec (x, y) = xy , donne 
d(Xf q X t 


rmn\ 

J 


(i) + M AT.-v, 


V 


mn\ 




V 

av (A T, X, 


mn > 
t y 


Piï (t) PT (t) ) d (B ij , B kl ) 


dxdy /t 

= A7 n a M (t) dt + X™ n fi pq (t) dBf + X™ n e^ kl (t) d (B ij , B kl ) t 
+x t v n (t) dt + AT7r dBf + AT77 (*) <* 

= X™ n dXf q + Xï q dX™ n + fd p k q (t) p% n (t) d (B ij ,B kl ) t . 


2. Provient des égalités suivantes : 


dB\ j dB kl = d(B ij ,B kl ) t , 
dtdBfi = dtd(B ij ,B kl ) t = 0. 


Remarque 1.2.3 Dans ce qui suit, on écrit (B l] ) t au lieu de (B l fi B l Q t . 

Définition 1.2.5 Comme dans le cas classique, on définit l’intégrale de 
Stratonovich o pour deux processus d’Itô matrices X et Y : 

X o dY = XdY + \dXdY et dX o Y = dXY + ^dXdY, 

où dXdY est le produit matriciel classique. 

La proposition suivante a lieu. 

Proposition 1.2.4 (voir [51],) Soient X et Y deux processus matriciels 
d’Itô. Alors on a 

( i) La formule d’intégration par parties : 

d (XY) = XdY + dXY + dXdY (1.14) 
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1 11 


d (XY) = dX 

O Y + X O dY 

(1.15) 

dXo(YZ) = 

(dX oY)o Z 

(1.16) 

(YZ)odX = 

(Y o Z) o dX 

(1.17) 

11 

h 

ST 

o 

= dY T o X T 

(1.18) 


et 


Y T désigne la transposée de la matrice Y. 

Preuve. Soient X, Y, Z des processus matrices d’Itô. 

(i) Pour tout i , j G 1 ,n, on a 

d (XY) ij = d ( X ik Y kj J =Y^ â C Xik Y kj ) 

= ( dX ik Y kj + X ik dY kj + dX ik dY kj ) 

k 

= ( dXY) ij + (. XdY) ij + (dXdYf , 


d’où 


d (XY) = XdY + dXY + dXdY. 

(ii) Par un calcul simple, on trouve 

dX o Y + X o d,Y = dXY + -dXdY + XdY + -dXdY 

2 2 

= dXY + XdY + dXdY = d (XY). 


On a 


dXo(YZ) = dX (YZ) + -dXd (YZ) 

= dX (YZ) + \dX (dYZ + YdZ). 

D’autre part 

(dX oY)oZ = (dXY + ^dXdY^j o Z 


dX (YZ) + ^dXdYZ + ^ fdXYdZ + ^dXdYdz) 
dX (YZ) + X -dX (dYZ + YdZ ), 
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1.2. G-MOUVEMENT BROWNIEN MATRICE ET G- CALCUL 

STOCHASTIQUE 

d’où (1.16). L’égalité (1.17) se démontre de la même manière. 

La formule (1.18) provient de la définition de l’intégrale de Stratonovich 
et du fait que 

(, dXdY) T = dY T dX T . 
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Chapitre 2 

G*—ED S des valeurs propres et 
des vecteurs propres des 
processus matrices 


L’objectif de ce chapitre est de donner les équations différentielles sto¬ 
chastiques des valeurs propres et des vecteurs propres des processus qui 
sont eux même solutions d’EDS gouvernées par un G—mouvement Brow¬ 
nien matrice. Pour contourner les difficultés liées au fait que les entrées d’un 
G—mouvement Brownien matrice ne sont pas forcément indépendantes, on 
supposera que leurs covariations quadratiques sont milles. 


2.1 Le modèle 

Dans le reste de cette thèse, on suppose que B satisfait l’hypothèse 
suivante : 

(A) Il existe un processus réel et croissant V tel que (B l i,B kl Y = SikSjib J t 
q.s. pour chaque i,j, k. I G 1, n, où ô uv est le symbole de Kronecker. 

On a alors a 2 t <b{< a 2 t où a := maxcr^ et a minerjj. Notons que 

i,j i,j - 

dans le cas classique, l’hypothèse (A) est satisfaite avec lP t — t pour tout 
t G [0, T\. 

On considère la G—équation différentielle stochastique générale définie 

par 


dX t = g (X t ) dB t h (X t ) + h (X t ) dB T t g (X t ) + o (X t ) dt + c (X t ) d (B) t (2.1) 
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2.2. G-ÉQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VALEURS PROPRES 


où les fonctions g, h, a, c : M —> M opèrent sur le spectre de X t , et X 0 G S n 
l’ensemble des matrices symétriques avec n valeurs propres différentes. La 
matrice variation quadratique d (B) est définie par d (B) := dBdB. Notons 
que grâce à l’hypothèse (A), d (B) est une matrice diagonale où d (B) lJ = 
SijdbL 

Dans ce qui suit, soit § n l’ensemble des matrices symétriques de dimen¬ 
sion n identifié à M"( n + 1 )/ 2 . Soit X t la solution de l’EDS (2.1) lorsqu’elle 
existe et soit H J X t H t — A t diag(\i(t )) la forme de diagonalisation 
de X t , où IIj est une matrice orthonormale. Notons que X G S n et que 
pour tout / : 1 -> 1 on a / (I) = Hf (A) H 7 . On suppose que X 0 G § n 
et que les valeurs propres de X t soient ordonnées de manière croissante : 

Ai (t) < A 2 (t) < ... < A n (t) jusqu’au premier temps de collision défini de 
la manière suivante : 

Définition 2.1.1 On définit le premier temps de collision des valeurs propres 
par 

t = inf {t : Xi (t) = A j ( t ) pour i fi j} . 


2.2 G —équations différentielles stochastiques 
des valeurs propres 

Maintenant on est en mesure d’énoncer le résultat principal de cette 
thèse. Notons que dans notre modèle, l’équation différentielle stochastique 
à étudier (2.1) se comporte comme dans le cas linéaire [4, 5,14,16, 37]. Les 
techniques utilisées sont inspirées par le cas linéaire, où le G— mouvement 
Brownien joue le rôle d’un mouvement Brownien classique. Pour cela on 
aura besoin de la proposition suivante : 

Proposition 2.2.1 [31] Soit X t G une solution de (2.1). Alors la ma¬ 
trice des valeurs propres A t est solution de l’équation différentielle stochas¬ 
tique suivante : 

dA = dN — dA o A + A o dA, (2-2) 

où dN = H 7 o dX o H et A est la matrice logarithme stochastique de H 
définie par 

dA = H 7 o dH. (2.3) 

Preuve. Comme dans [4,16], on considère la décomposition 

A = H t XH, 
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CHAPITRE 2. G-EDS DES VALEURS PROPRES ET DES 
VECTEURS PROPRES DES PROCESSUS MATRICES 


d’où 

dA = d ( H t XH ) . (2.4) 

D’après la formule (2.4) et les notations d’Itô Stratonovich, on peut écrire 

dA = dH T o (XH) + H t o d (XH) 

= (d,H T o X) o H + H t o (dX oH + X o dH) 

= ( dH T o X) o H + H t o dX o H + H t o X o dH , 

d’où en posant dN = H T o dX o H, 

dA = ( dH T o HAH t ) oH + dN + H T o HAH T o dH. 

On obtient en utilisant les formules (1.16) et (1.17) : 

dA = dN+ dH T o HA +H T HAH T odH 
= dN+ dH T oHA + AH T odH 
= dN + dH T oHoA + AoH t odH 
= dN — dA o A + A o dA. 


Théorème 2.2.1 [51] Soit X t une solution de l’équation (2.1) telle que 
X 0 £ S n . Alors il existe un G—mouvement Brownien réel w l tel que = 
8ijV. Pour tout i,j £ 1 ,n et pour t < r le processus des valeurs propres A f 
est solution du système suivant : 

d\i = 2 g (Ai) h (Ai) Y,H ki dw k + a (Ai) dt + dV ü (2.5) 

k 

OÙ 

dV ij = 8 ijC (Ai) £ [H ki ) 2 db k + dR ij (2.6) 

k 

et 

dR ij = [<W(A fc )ù(A fc )MA,-) (2-7) 

k^j N k 

+g 2 (\ k )h(\ i )h(\ j )'EH li H lj db l 

i 

+6, j g 2 (\,)h 2 (\ t )Y.(H ,k Ÿdb l . 

I 
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2.2. G-ÉQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VALEURS PROPRES 


Preuve. Par souci de simplicité, on utilise la notation d r ' J au lieu de 
(X t y j pour d' = g, g 2 , h, a et c. D’après la proposition précédente on 
d\ = dN — (IA o A + A o (IA. Observons maintenant que le processus 
o dA — dA o A est nul sur la diagonale. Par conséquent dX, = dN n et on 
a 0 = dN l i + (Ai — A j) dA si i y j et donc 

dA zj = -——— dN zj pour i y j (2.8) 

Xi 

On a d’après l’équation (2.1) 


dX ij = £ g ip dB pq h qj + E h ip dB qp g qj + a ij dt + E c ip d < B pq , B qj ) f , 


p,q 


p,q 


p,q 


d’où 

dxi j dx£ m 


£ g ip dB pq h qû + £ h ip dB qp g qj 
p,q p,q y 

* f £ g kp 'dB p q ' h q m + £ h kp 'dB q ' p ' g q '‘ 
\p',q' p,q 

Y g ip dB pq h qj Y g kp 'dB p ' q 'h q ' m 


p,q 


p ,q 


Y g ip dB pq h qj Y h kpl dB q,p 'g q ' m 


p,q 


+ Y h ip dB qp g qj Y g kp 'dB p,q 'h q ' m 

p,q p',q' 

+ Y h ip dB qp g qj Y h kp 'dB q ' p 'g q ' m 

p,q p',q' 

= E g ip h q ig kp h qm 8 pp ,8 qq ,db q t + E g ip h qj h kq g pm 8 pq ,8 qpl db q 
p,q,p',q' p,q,p',q' 

+ E h ip g qj g kq h pm 8 qp ' 8 pq 'db p + E h ip g qq h kp g qm 8 qql 8 ppl dlf t . 

p,q p,q 

Il résulte que 

dx l t j dx km = Yy\g iVhqj g kPhqm + g iphqj g pmhkq }db q t 

p,q 

+ E [h ip g qj g kq h pm + h ip g qj h kp g qm ]dbï 


p,q 


= E {g 2 ) ik h l ih lm +(g 2 y m h l jh kl 


+ (g 2 ) kj h il h! m + (g 2 ) jm h il h kl 


db\ 
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VECTEURS PROPRES DES PROCESSUS MATRICES 


Finalement, on obtient 


dX\ j dX km = £ [ G g' 2 y k h lj h lm + (g 2 ) kj h li h lm 

i 

+ (g 2 Y m h l ^h lk +(g 2 Y m h li h lk 1 db\ 


Il est facile de vérifier que 


d,N = H T dXH + -H T dXdH + - dH T dXH, 
2 2 


( 2 . 10 ) 


il résulte alors que la partie G—martingale de dN est égale à la partie 
G—martingale de II r dXII. donnée par : 


dN ij dN km = (H T dXH) ij (H T dXH) km 

= H pi dX pq H qj ^ H p ' k dX p ' q ' H q ' r 


(2.il) 


La formule (2.11) écrite avec (p,q,p',q') au lieu de ( i,j,k,m ), devient 
dN ij dN km = E E H pi H qj H p ' k H q ' m (g 2 ) PP ' h lq h w + (g 2 ) p ' q h lp h lq ' 

l p,q,p' ,q' 


(tg 2 ) pq ' h lq h lp '+ (g 2 ) qq 'h lp h lpl } db 


( 2 . 1 . 2 ) 


= E E H pi (g 2 ) pp H pk J2H qj h lq h lq H qm db 1 

l p,p’ q,q' 


+E E H ca (g 2 ) P q H p ' k E H pt h lp h lq 'H q ' m db 1 

l q,p' q',p 


+E E HPl (f/ 2 )” Rq ' m E H<ü h lq h lp 'H p ' k db 1 

l p,q' q,p' 


+E {g 2 ) qq 'H q ' m E HPi h lp h lp 'H p ' k db 1 

l q,q' p,p' 
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2.2. G-ÉQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VALEURS PROPRES 


Observons maintenant que 

E H- {gY H*' k = (H T fH) ik = g ' 2 (A)* = 8 ik g 2 (X k ), 

p,p' 

d’où 

(/) = WMeWe^^W"* 

q,q' V i J 

= S ik g 2 (À fc ) E (h (X) d (B) h (X))™' H«' m 

q,q' 

= 5 ik g 2 (A fc ) (J H T h (X) d (B) h (X) H) jm , 
et de façon similaire, on a 

(II) = 8j k g 2 (X k ) ( H T h (X) d (B) h (X) H) im , 

(III) = S im g 2 (A m ) ( H T h (X) d (B) h (X) H) jk , 

(IV) = 8 jm g 2 (X m ) {H T h (X) d (B) h (X) H) ik . 

Il résulte que 

dN ij dN km = ô ik g 2 (Àfc) [H T h (X) d (B) h (X) H) jm 
+S jk g 2 (Afc) (H T h (X) d (B) h (X) H) im 
+S im g 2 (X m ) ( H T h (X) d (B) h (X) H) kj 
+8 jm g 2 (X m ) ( H T h (X) d (B) h (X) H) ik , 

et en utilisant le fait que h (A") H = H h (A), on obtient 

dN ij dN km = S ik g 2 (Afc) (h (A) H T d (B) Hh (A)) jm 
+8 jk g 2 (Afc) (h (A) H T d (B) Hh (A)) im 
+S im g 2 (X m ) (h (A) H T d (B) Hh (A)) kj 
+S jm9 2 (X m ) {h (A) H T d (B) Hh (A)) lk . 

Comme (H T d(B) H)‘ 2 = E H h H l ^db\ alors 

i 

(h (A) H T d {B) Hh (A)) ij = h (Xi) h (Xj) E H li H lj db l , 


(2.13) 
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d’où d’après la formule (2.12) 

dN^dN km = 8 ik g 2 (À*;) h (Xj) h (À m ) E H lj H lm db l (2.14) 

i 

+S jk g 2 (Afe) h (A,) h (À m ) E # '"‘dô 1 

i 

+S im g 2 (A m ) h (A fc ) h (Aj) E # y d&* 

i 

+S jm g 2 (A m ) h (Ai) h (X k ) E # “i^dô 1 , 

i 

et donc, 

dN u dN jj = 4:8ijg 2 (Ai) ù 2 (Ai) E (H li ) 2 db l . (2.15) 

i 

Il s’ensuit qu’il existe un G— mouvement Brownien w l satisfaisant (u;\ uù) = 
8ijV , tel que la partie G—martingale de c/A* soit égale à 

2g (Ai) h (Ai) E H ki dw k (2.16) 

k 

Passons maintenant à la partie à variation finie dF de dlV. On a 

dF t = H T aHdt = a (A t ) dt. 

On en déduit que F est diagonale et 

dFf = a(\i(t))dt (2.17) 

D’autre part, la partie intégrale dV par rapport à dlf de d,N : 

dV = H T cd (B) H + ^ (dH T dXH + H T dXdH) 

= H T cHH T d (B) H + ^ (( dH T H ) (. H T dXH ) + ( H T dXH ) ( H T dH )) 
= c (A) (F) F + ^ (dNdA + (dNdAf^j . 

Notons que (c (A) H T d (B) H) 1J = ôijC (Ai) E (H kl ) 2 db k . Il s’ensuit que, si 

k 

(,INdAf = E dN ik dA kj = E 1 dN ik dN kj . 

k k^j A j X k 
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2.2. G-ÉQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VALEURS PROPRES 


On a d’après la formule (2.14) avec (i, k, k,j ) au lieu de (i,j, k , m), 

1 


(dNdAy j = y 


k^j Xj — A k i 


<W (A fc ) h (X k ) h (Xj) YH lk H lj db l 


+g 2 (X k )h(X i )h(X j )EH li H lj db l 


+S ijg 2 (Xj) h 2 (Afc) YH lk H lk db l 


= E 


+4;/ (Aj) h (Xi) h (X k ) YH li H lk db l 

i 

' (Afc) h (X k ) h (Xj) YH lk H l >db l 

i 


k^j Xj X k |_ 


V(Afc)/i(A i )/ l (A i )E^ y ^ 

i 


et 


(dNdAf = Y 


k^i Xi X k |_ 


g 2 (Afc) h 2 (Aj) Y (H li ) 2 db l 


+g 2 (X l )h 2 (Xk)Y{H lk Ÿdb l 


ce qui implique que, pour i y j 


dV ij 



j_ 

—Afc 


et que 


(W (Afc) h (Xk) h (Xj) + g 2 (Afc) h (Xi) h (A,)) YH h H l >db l 

i 

(2.18) 


dE* = c (Ai) E 2 db k + E 


k^i 




+9 2 (A,)ft 2 (A t )£(ff“f Æ/ 


g 2 (\ k )h 2 (\,)j:(H l -) 2 dl , 1 

l 

(2.19) 


Par conséquent, d.Xi = 2g (Ai) h (Ai) E H kl dw k + a (Ai) dt + dV u , d’où le 

k 

résultat désiré. ■ 

Notons que dans [34, 35], les auteurs traitent le cas des mouvements 
Browniens avec des méthodes analogues. 
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CHAPITRE 2. G-EDS DES VALEURS PROPRES ET DES 
VECTEURS PROPRES DES PROCESSUS MATRICES 


2.3 G —équations différentielles stochastiques 
des vecteurs propres 


Pour expliciter les G—équations différentielles stochastiques des vecteurs 
propres, on aura besoin de la proposition suivante : 

Proposition 2.3.1 [31] La matrice H t satisfait l’équation suivante : 

dH = H o dA = HdA + -HdAdA. (2.20) 

2 


Preuve. On a 


dA = H T odH 

= H T dH + ^ d,H T d,H. 

En multipliant cette égalité à gauche par H et en ajoutant \dHdA des deux 
côtés, on obtient 


HdA + ^dHdA = dH + ^HdH T dH + ^ dHdA. 


D’une part, on a 


- ( HdH T dH + dHdA) = - (HdH T dH + dH (H T dH + -dH T dH 


2 V 

= ^ ( HdH T + dHH T ) dH, 

d’autre part, on a 

0 = dl = d ( HH t ) = HdH T + dHH T + dHdH T , 

d’où 

Hd,H T + d,HH T = -dHd,H T , 

et 

^ ( HdH T dH + dHdA) = ~dHd,H T dH = 0, 

et en combinant ces égalités, on obtient 


dH = H odA = HdA + -HdAdA. 
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2.3. G-EQUATIONS DIFFÉRENTIELLES STOCHASTIQUES DES 

VECTEURS PROPRES 


Théorème 2.3.1 (voir [51] ) La matrice des vecteurs propres H t satisfait 
l’équation différentielle stochastique suivante 


dHij = Ettt- 

k^j Aj Ak 


g (A k ) h (Xj) ( d/3H) jj + g (A,) h (A fc ) (d^H) h 


+dV kj } - - £ H ik dQ kj 

2 k 

où dV k i est donné par l’équation (2.6), 


( 2 . 21 ) 


dQ k> = E TT-TTTT-TT hiS 1 (AJ h 2 (A,) E (H*)-‘d.v 

l^k,l^j l Al — A k ) [Al — Aj) l p 

+g 2 (Ai) h (A fc ) h (A,) ZH pk H^dlf\ , 


et (/3 ÎJ ) est un G—mouvement Brownien matriciel satisfait l’hypothèse (A). 
Preuve. On a d’après la proposition (2.3.1) 


d.H = HdA + -HdAdA, 


dH ij = H ll dA lj + -^Th ü ( dAdA) lj 

i i 

D’une part, on a d’après l’équation (2.8) : 

(dAdAf = JfdA ik dA kj 

k 

= E U- A HA ,P Ni “AN ki 

k+i,k^j [Ak — Ai) [Aj — Ak) 


( 2 . 22 ) 


= E TT—rrn—n %/ 2 (A,)fc 2 (A t )E(»'T<«>'+ 

k^i,k^j \A k — Ai) t Aj — A k ) L l 


g 2 (Afc) h (Aj) h (Xj) H li H lj db l 


D’autre part, à l’aide de l’équation (2.14), on déduit que si i j, 


dN ij dN ij = g 2 (A*) h 2 (A j) E (H lj ) db l + g 2 (Xj) h 2 (A f) £ (H li ) db l . 

i i 
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CHAPITRE 2. G-EDS DES VALEURS PROPRES ET DES 
VECTEURS PROPRES DES PROCESSUS MATRICES 


Il s’ensuit que la partie G—martingale de dN' :) est 


9 (A*) h (A,) E H lj d/3 jl + g (Xj) h (Xi) E H li d(3 ü 

i i 

g (Xi) h (Xj) (d(3H) jj + g (A,) h (A,) (d(3H) ü ] , 


où (3 := (T 3 ) est un G—mouvement Brownien matriciel satisfait l’hypo¬ 
thèse (A) et que si i ^ j 


dA ij 


1 

Xj Xi 

1 

Xj A i 


dN ij 

g (Xi) h (A,) (d(3H) jj + g (A,) h (Xi) (df3H) ü 


(2.23) 

+ dV ij 


La formule (2.21) découle de (2.20), (2.22) et (2.23). ■ 


2.4 Temps de collision 

Le Théorème de Lévy ci-dessous nous donne une manière facile d’af¬ 
firmer quand est-ce-qu’une martingale locale continue est un mouvement 
Brownien. 

Théorème 2.4.1 (Théorème de Lévy voir [47] J 

Soit W une martingale locale continue de dimension n x n, nulle en 0 
telle que 

/ppù W^\ _ / si i k et j / 

' ’ 3t \ 0, sinon 

Alors W est un mouvement Brownien. 

Preuve. Il suffit de montrer que X := vec(W) est un mouvement Brow¬ 
nien n 2 —dimensionnel. Soient a, b G {l,...,n 2 }. Alors, il existe i,j,k,l G 
{1,..., n} tel que a — n (j — 1) + i et b — n (l — 1) + k. Clairement, i = k et 
j = l impliquent a = b, et inversement. Ainsi, 

(X a ,x b ) = (W i fW kl )=tl {i=k} l {j=l} =tl {a=b} . 

On conclut d’après le théorème 40, chapitre 2 Protter [48] que X est un 
mouvement Brownien n 2 — dimensionnel, u 
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3.1. THÉORÈME DE YAMADA-WATANABE 

MULTIDIMENSIONNEL 


( i ) Il existe une fonction p strictement croissante sur M + , nulle en 0 et telle 

que - p ^ L 1 (M + ) et \a ( x ) — cr (y)\ 2 < p (|æ — y\) pour tout 

(ii) b est une fonction continue Lipschitzinne. 

Alors l’EDS 

dX t = a (X t ) dB t + b (X t ) dt (3.1) 

admet une solution forte unique. 

P. Graczyk et J. Malecki sont les premiers auteurs qui ont proposé une 
version multidimensionnelle de ce théorème, qui a été la clé de l’étude de 
l'existence et l’unicité de l’EDS de type (2.1). On aura besoin du lemme 
suivant 

Lemme 3.1.1 [21,30] (Lemme de Gronwall) Soient p une fonction conti¬ 
nue non négative sur un segment [a, b] et vérifiant l’inégalité : 

t 

Vt G [a, b] ; p (t) < C + A f p (s) ds. 


Alors 

Vi G [a, b] \p{t) < Ce At . 

Théorème 3.1.2 [16] Soient p,q, r G N et ai : W —■> M, i = 1 ,p et 
Ck,dj : M p+r —» R, k — p + 1, ...,p + q ; j = p + 1, ...,p + r des fonc¬ 
tions à valeurs réelles bornées et continues, satisfaisant les conditions de 
Lipschitz suivantes : 

1- k (l/i) - ai ( 2 / 2)1 < A ||i/i - 2 / 2 II — 

2. \c k (y 1 ,z 1 ) - c k (y 2 , z 2 )\ < A^y^zf) - ( 2 / 2 , ^ 2 ) || ,k = p + l,p + q; 

3. \dj ( 2 / 1 , ^ 1 ) - dj (y 2 ,z 2 )\ < AWfy^zÉ) - ( 2 / 2 , z 2 ) ||, j =p+l,p + r, 

pour tout yi,y 2 G M p et z\,z 2 G M r . De plus, soit cq : M —> M,i = l,p 
un ensemble des fonctions boréliennes bornées telles que 

Wi(x) ~(Ji(y)\ 2 <Pi(\x-y\), x,yeR 

où p i : (0, 00 ) —> (0, cxd) sont des fonctions boréliennes telles que : 

Jpf 1 ( x ) dx = 00 . 

0 + 
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CHAPITRE 3. EXISTENCE ET UNICITÉ DE LA SOLUTION 
D’UNE G-ÉQUATION DIFFÉRENTIELLE STOCHASTIQUE 

MATRICIELLE 


Alors le système des équations différentielles stochastiques suivant 

dY i = ai (Y 1 ) dB i + a t (y) dt; i = ïfp (3.2) 

p+q 

d,Z j = Ck (Y, Z) dB k + aj (Y, Z)dt\j — p + l,p + r (3.3) 

k—p+l 

admet une solution forte unique. 

Signification de p, q et r : 
p : nombre de coordonnées de Y, 

q : nombre de mouvements Browniens dans les EDS pour dZ, 
r : nombre de coordonnées de Z. 

Preuve. La preuve est basée sur l’approche présentée dans Revuz,Yor [49], 
en particulier sur les résultats de Le Gall [24], 

Soit {Y, Z) et (y, Z^j deux solutions de (3.2) et (3.3) définies avec le 
même mouvement Brownien B = {Bf) i<p+q telles que : 

Y 0 = Ÿ 0 et Z 0 = Z 0 p.s.. 


On a pour tout i — 1, 


t 

Y i -Ÿ i = J (Yj) - a, : (î?)) dB\ + j ( a, (Y s ) - a t [Y,) )ds (3.4) 
o o 


ce qui implique que 


d{Y l - r\Y‘ - Y‘) = (<7, (Yi) - a, (pj J (a, (Y<) - a, (p j j dt 

(3.5) 

t 

f y s>ÿ£ 

et en multipliant par / , on obtient 

I Pi ( *s ~ ) 


l {Yi>Yj} 


Pi [Y*-Y* 


d(Y i -Y i ,Y i -Y 1 ) = 


o 

< t, 


( g < ( y ;)- g < ( î ?)) 2 
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3.1. THÉORÈME DE YAMADA-WATANABE 
_ MULTIDIMENSIONNEL 

d’où d’après le lemme 3.3 de [49] p.389, le temps local de Y 1 — Y 1 en 0, 
I}] (Y 1 — v'j est identiquement nul. Par conséquent, par la formule de 
Tanaka on a 

t 

P - p| = y sgn (y; - p) d (y; - p) + (y - P) 

0 

t 

= [ sgn {y; - P) d (y; - p) 

0 

t 

= f sgn (y; - P) ((T. (P) - (T. (P) ) dBl 
0 

t 

+ J sgn (Y; - y;) («* OY - a* (Y)) 

o 

Comme a, est bornnée, alors 

t 

Y ~ Y| - j sgn (Y - Y) (y) - a, (Y)) ds 

o 

est une martingale nulle en 0, d’où d’après les conditions de Lipschitz sa¬ 
tisfaites par les a*, 

t 

e Y - Y - A J E Ys ~ Y ds - 

o 

En sommant ces dernières inégalités par rapport à i. on obtient 

t 

£|y-Y|| <cJ e\\y s -Ÿ 8 jds, 

o 

d’où d’après le lemme de Gronwall Y t = Y t pour tout t > 0 p.s.. 
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CHAPITRE 3. EXISTENCE ET UNICITÉ DE LA SOLUTION 
D’UNE G-ÉQUATION DIFFÉRENTIELLE STOCHASTIQUE 

MATRICIELLE 


Il résulte, d’après les propriétés de l’intégrale d’Itô et de l’inégalité de 
Schwartz, que pour tout t G [0; T] 

2 

«J -r, / <■ \ 

12 


E 


zi - zi 


< 


P+q ! r. 

C J2 E [ {c k (Y s ,Z s )-c k (Y s ,Z s ^dB k s 

k = P + î '^ 


+CE I [a, (Y s , Z s ) - dj Y s , ZA\ds 


p+q , 

< G Y. E J ( Ck (Y ^ Zs 


) - c k [Y S1 Z s ) ) ds 


k=p+l 0 
t 


ACTE / cij (Y a , Z a ) - aj Y s , Z s )) ds 


< CA 2 (q + T)E 


Y s -Y s 


+ 


Z s — Z s 


ds. 


Ainsi, d’après l’argument utilisé précédemment on a Y = Y et on obtient 


E 


z t -z t 


< CA 2 (q + T) r E 


Z s — Z s 


ds. 


On conclut d’après le lemme de Gronwall que Z t = Z t p.s. pour tout t G 

[o,n- ■ 


3.2 Cas des G— équations différentielles sto¬ 
chastiques 

Dans ce paragraphe, on utilise les hypothèses de croissance linéaire et 
de Lipschitz sur les fonctions g, h, a, c comme fonction de matrices (voir 
l’appendice) et non de variables réelles. 

Hypothèses : On suppose que g, h, a, c : M — > M, satisfont les condi¬ 
tions suivantes : 

(i) Condition Lipschitz : pour tout X , Y G M nxn 

| J(X) -J (Y) | 2 < A\X - Y \ 2 , 
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3.2. CAS DES G-ÉQUATIONS DIFFÉRENTIELLES 

STOCHASTIQUES 


où J (X) = g (. X) tk h (. Xy l , a (X) et c (X) respectivement et A est un 
constante positive. 

(ii) Condition de croissance linéaire : pour tout X G M nxn 

|J(A')| 2 < A'(l + |Af), 

où J (X) = g (X) lk h (Xy l , a (X) et c (X) respectivement et K est 
un constante positive. 

Théorème 3.2.1 [51] Sous les hypothèses (i) et ( ü ), on a l’unicité en 
trajectoires de la solustion de l’équation différentielle stochastique (2.1). 


Preuve. On a 

dX ij = E(g ( x Ÿ k h ( x ) lj + 9 C x Ÿ j h ( X ) Ü ) dB kl +a (.X) ij dt+c (X) ij d (.B jj ) 

k,i ' ' 

(3.6) 

Soient X (xf) ,k = 1,2 deux solutions de l’EDS (2.1) avec les conditions 
initiales Xk = (x'f ). Par souci de simplicité, on utilise la notation 

jijH pQ , = g ç x yk h ç x yj + g ^ X jkj h (xf 


pour X G M nxn . Alors on a pour tout u < t 

xÿ (*o - xÿ (* 2 ) = \xf-xi j \ + 


(■ J ijkl (x 8 (xf) - r kl (X s (x 2 ))) dB 

k,l 


kl 


+J (a (■ X s (xi)) ÎJ - a (X s (x 2 )Y J j ds 
o 

u 

+ [ (c (X s ( Xl )) ij - c (X 8 (x 2 )f) d (B ”) S , 


d’où 


E Xif (x,) - Xij (x 2 ) <c\ I X? - X? + (!') + (II') + (III') , 


ou 


V) = E E 


k.l 


(II') = E 


(J ljkl (X s (xf) - J ijkl (X s (x 2 ))) dB 


a (X 8 (xi)y j - a (X 8 (x 2 )Y 3 ) ds 


kl 
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CHAPITRE 3. EXISTENCE ET UNICITÉ DE LA SOLUTION 
D’UNE G-ÉQUATION DIFFÉRENTIELLE STOCHASTIQUE 

MATRICIELLE 


(///') = E ( (c (X s ( Xl )f - c (X s (x 2 )ÿ j ) d (B 33 ) 


Pour (/'), on applique l’inégalité BDG (1.9) avec p — 2 et on a 


(/') < E I supV / {,T 3kl {X s (xQ) - J ijkl (X s (x 2 ))) d,B kl 

y k,i q 

< e(j2 sup f ] {J ijkl (X s (xQ) - .T 3kl (X s (x 2 ))) dB k s 


W V o 


< l E (\ jijkl ( Xs ^ - jijH (^fds) > 


d’où d’après la condition de Lipschitz 


(/') < C E £(|X s (z 1 )-X s (z 2 )| 2 ds) 


< CEI sup |X U (xq) - X u (x 2 )| ds 

J \u<s 

0 


Pour (//') 


(//') < £ / ds 


< E ( \X S (xQ 13 - X s (x 2 ) ?J ) ds 


< C E [sup \X U (xi) - X u (x 2 )| ds . 
J \u<S J 

0 
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3.2. CAS DES G-ÉQUATIONS DIFFÉRENTIELLES 

STOCHASTIQUES 


On a d’après l’inégalité BDG (1.8) pour p — 2 


(IIT) < et I e[\(c(x s (; Xl )f - c (X s (x 2 )) i 


ds 


,0 


< et | JE 

,0 


x s ( Xl y j - x s ( X2 y 


d,s 


< CT ( El sup \X V (xQ - (x 2 )| 2 ds . 


u<s 


On obtient finalement en sommant (/'), (IT) et (III') 


E sup \X U (xi) - X u {x 2 )\ 2 ) < C (T, n) (|xi - x 2 \ 2 


U<t 


+ / E ( sup \X U (xi) - X u (x 2 )\ ) ds | , 

J \u<S 

0 


d’où d’après le lemme de Gronwall, 

0 < E (sup \X U (xi) - X u (x 2 )| 2 ^) < C (T, n) \ Xl - x 2 \ 2 e c(T ’ n)t . 

\u<t J 

En particulier, si X\ = x 2 on obtient l’unicité en trajectoires de X t . ■ 

Théorème 3.2.2 [51] Sous les hypothèses (i) et (ii) VEDS (2.1) admet une 
solution unique. 

Preuve. Pour l’existence de la solution de l’EDS (2.1), on considère la suite 
de Picard m X = ( m X l3 ) m&i définie par 

°X l t j — x ij G M, 0 < t < T 


et pour wéN 

t t t 

171+1 X^ = x ij +J2 J J ijH ( m Xs) dBf+ J a ( m X s ) ij ds+ J c Ç m X s ) ij d < B jj ) g 
k ’ 1 0 0 0 

(3.7) 
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CHAPITRE 3. EXISTENCE ET UNICITÉ DE LA SOLUTION 
D’UNE G-ÉQUATION DIFFÉRENTIELLE STOCHASTIQUE 

MATRICIELLE 


On a 


E ( \ m+1 X \ j | 2 ) < C'I \x ij \ 2 + 


k,l 


J ijkl ( m X s ) dB 


+E 


+E 


a ( m X s ) ij ds 


c ( m X s ) ij d 


Par un raisonnement similaire que celui utilisé dans le théorème précédent, 
on obtient d’après l’hypothèse de croissance linéaire : 


E ( \ m+1 Xi j 2 ) < C (T, n) I x ij z + / ( 1 + E( m Xl ] V ))ds , 


et en sommant par rapport à % et j. on obtient 


E 


(r%| 2 ) <C'{T,n) l\x\ 2 + T+ f E (\ m X s \ 2 ) ds . 


Il résulte que, d’après le lemme de Gronwall, 

E ( \ m+1 X t \ 2 ) < C (T, n) (\x\ 2 + T) e c ' (T ’ n)t . 

Maintenant, on montre que (' m X ) est une séquence de Cauchy dans L 2 G (Dr ) 
On a 


k+1 XÎ j - k XÎ j = ^ (J ijkl ( k X s )-J ijkl ( k ~ 1 X s ))dB k 

O o 

t 

+f (a ( k X s ) ij - a (^ 1 X S ) !3 ) ds 

0 

t 

+ [ (c( k X,f-c( i - I X 
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3 . 2 . CAS DES G-ÉQUATIONS DIFFÉRENTIELLES 

STOCHASTIQUES 


Par un argument similaire à celui utilisé dans la preuve de l’unicité, on 
obtient que 


E ( k+ 1 X t - k X t 2 ) < C (T, n) ( Ie( k X s - 2 ) ds I ( 3 . 8 ) 


< c(T, n y 



E ( l^x 


i 2 - k ~ 2 X t2 2 ) dudt 2 , 


et en répétant cette procédure on obtient 

t ti ifc 



E ( k+ 1 X t - k X t \ 2 ) < C (T, n ) k+1 / / • • • x X tk - x 2 ) dt^h ■ ■ ■ dt 



0 0 0 


D’autre part, on a 


EiVX^-xl 2 ] = E 



tk tk 

J2 J J ijpq (x) dB p s q + J a (x) ij ds 


p,g o 


^k 


+ c(x) i ’d{B») < 


! \ 


) 


J ijpq B pg + a ( x )ÿ t k + c ( x ÿi (Bü} 


tk 


i.j \ p,q 


< 


K (n) Y, 


1,3 


Y{J iir 'MŸE(B”Ÿ+(a(xf) T 


p,g 


+ (cW«) 2 e((s*>,J ; 

< K(n,x,â,T), 


d’où 


\ k+1 X t 


*X,\\l = E (|‘ +1 V - ‘V,| 2 ) < K(n,x,c,T) E 


fe+1 
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CHAPITRE 3. EXISTENCE ET UNICITÉ DE LA SOLUTION 
D’UNE G-ÉQUATION DIFFÉRENTIELLE STOCHASTIQUE 

MATRICIELLE 


et par conséquent pour tout p, m G N 

m-\-p— 1 

,2 < 


\\ m+p x t 


171 X 


t O ^ 


k=m 


< v ir% - “x.i < ^ 

k=m 

oo 

< yjK (n,x,W,T) 


k X t 


(C” ( T, n ) T) 


k +1 


k=m 


(k + 1)! 


On conclut que ( m X t ) est une suite de Cauchy. 

Soit maintenant X, la limite de m X t dans fQ /d et montrons que X 4 
est bien la solution de l’équation (2.1). Il suffit alors de prouver que 


lim E 

m—>oo 



Jijpq (X s )) dB Pq 


lim E 

m— kx) 




ds 


= 0, 


= 0 


et 


lim E 

m— kx) 



c ( m x s y j - c (x s y j )d{B”) =o, 


La première et la troisième égalité sont assurées par les conditions de Lip- 
schitz et l’inégalité BDG. Pour la seconde égalité, on a d’après l’inégalité 
de Hôlder et la condition de Lipschitz : 



a (X.f 



< TA 


l X q -XMds. 


On conclut en appliquant les deux membres de cette inégalité E et le fait 
que m X converge vers X. m 
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Chapitre 4 
Appendice 

4.1 Algèbre matricielle 

Dans cette section, on résume les définitions et les techniques préalables 

qui sont nécessaires pour les matrices aléatoires. 

Définition 4.1.1 

(i) § n := le sous-espace linéaire de toutes les matrices symétriques de 

M n (R). 

(n) §+ ( resp . §“) = l’ensemble de toutes les matrices symétriques définies 
positives (négatives) de A4, n (M). 

(iii) §+ = la fermeture de dans M. n (R), c’est-à-dire l’ensemble de 
toutes les matrices symétriques semi-définies positives de A4 n (M). 

Le théorème suivant décrit les matrices définies positives. 

Théorème 4.1.1 (Matrices définies positives) 

(i) A G S+ si et seulement si x T Ax >0VïGR n :i^0. 

( ii) A G §+ si et seulement si x T Ax > 0 V i G R n : ||a;|| = 1. 

{iii) A G S+ si et seulement si A est diagonalisable par une matrice or¬ 
thogonale avec des valeurs propres positives, i.e. il existe une mutrice 
orthogonale U G A4 n (K), UU 1 = I n , telle que A = UDU T avec une 
matrice diagonale D = diag \ n ), où Xi > 0, i = 1 sont 

les valeurs propres positives de A. 

( iv) Si A G §+ alors d” 1 G §+. 

(v) A T A G §+ pour toute matrice A inversible. 
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(■ vi) A G §+ si et seulement si A est diagonalisable par une matrice ortho¬ 
gonale avec des valeurs propres non négatives. 

{vii) M t M G §+ pour tout M G M n (M). 

Preuve, (voir Muirhead [32] ou Fischer [11]). ■ 

Définition 4.1.2 Soit S G A4 n (M) . On définit l’opérateur différentiel Dÿ = 
(q§ïj)‘ J pour toutes les fonctions différentiables f : A4 n (M) —> M, la matrice 
de toutes les dérivées partielles Dijf (S) de f (S). 

Règles de calcul pour les déterminants 

Pour tout A,B e .M n (R), S G GL(n) l’ensemble des matrices inver¬ 
sibles et H t : M —> GL (■ n ) différentiable, alors 

Lemme 4.1.1 

(i) det (AB) = det (A) det (B) et det (aA) = a n det (A) Va G M. 

(ii) Si A G GL (n) ou B G GL (n) alors det (I n + AB) = det (I n + B A). 

(iii) | det (H t ) = det (H t ) tr (Hfi l f t H t ) . 

(iv) D (det (S)) = det (S) (S. 

(v) det (A) est le produit des valeurs propres de A. 

En outre, si S est symétrique, alors 

(vi) D (det (S)) = det (S) SS 

(S^f (S~ 1 ) ij - (S- 1 )^ (S^f 

Où (5 ,_1 ) y désigne la i, j—ième entrée de F -1 . 

Preuve, (voir Fisher[ll] et [47,49]). ■ 

Lemme 4.1.2 (Règles de calcul pour la trace). 

Pour tout A, B, S G A4 n (M) , 

(i) tr (AB) = tr (BA) et tr (aA + B) = atr (A) + tr (B) Va G M. 

(ii) tr (A) est la somme des valeurs propres de A. 

Preuve, (voir Fischer [11]). ■ 

La définition suivante nous donne une relation entre les vecteurs et les 
matrices. L’idée est la suivante : supposons que nous voulions transférer 
un théorème qui tient pour des processus stochastiques multivariés à un 
qui tient pour des processus stochastiques matriciels S. Ensuite, on peut 
appliquer le théorème à vec (S) aux processus multivariés résultants pour 
obtenir le théorème dans une version matricielle. 


(vii) QS if ds u (det (S)) = det (S) 
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Définition 4.1.3 Soit A G A4 m ,n (M) avec les colonnes ai G = 

1, on définit la fonction vec : Ai m ^ n (M) —■> W' nn par 


vec (A) 


/ ai \ 

y J 


Parfois, nous considérerons aussi vec (A) comme un élément de A4 mnjl (M) 


Remarque 4.1.1 


vec 


P T ) = 


/st\ 


V < / 


où Oj G M n , j = 1, ..., m désigne le vecteur colonne de A. 


Lemme 4.1.3 (voir[ 18] Théorèm {1.2.22)) Pour A, B G M m , n (R) il est 
vrais que tr ( A T B ) = vec{A) T vec{B). 


4.2 Matrices aléatoires 

Dans cette thèse, on a utilisé les définitions sur les processus matrices 
suivantes : 

Définition 4.2.1 Une fonction mesurable X : R + x D —> M. mn (R), 
{t, w) —> A" {t, w) = X t (w) est appelée un processus stochastique (variable 
matricielle) si X ( t,w ) est une matrice aléatoire pour tout t G M + . 

Définition 4.2.2 Un processus stochastique matriciel X est appelé martin¬ 
gale locale, si chaque entrée de X est une martingale locale, i.e. s’il existe 
une suite croissante de temps d’arrêt, strictement monotones {T n ) n£N , T n 
oo, tel que {X tr , Tn ) l: j est. une martingale pour tout, i, j. 

Définition 4.2.3 Un mouvement. Brownien matriciel B dans M. n p (M) est. 
une matrice composée de mouvements Browniens réels indépendants, i.e. 
B = {B l] ) i 3 où B' 3 sont, des mouvements Browniens indépendants et. uni¬ 
dimensionnels, 1 < i < n, 1 < j < p. 

Définition 4.2.4 Un processus stochastique matriciel X est. appelé semi- 
martingale si X peut être décomposé en X = X 0 + M + A où M est. une 
martingale locale et A un processus adapté à variation finie. 
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4.3 Fonction de matrices 

Soit des entiers strictement positifs et n leur somme. Soit 

Ai,A r des nombres complexes distincts. On pose 

r 

p 0 (x) = H(x-x l r 

i— 1 

Soit A une matrice dont le polynôme minimal est P 0 , et r un endomor¬ 
phisme dont la matrice est A dans une base donnée. 

Dans la suite, on pose £p 0 = C [X] / Pq [X] et A4 a = {P (A) \P G C [A"]}. 

Définition 4.3.1 Soient A et B deux matrices de C n,p , on dit que A est 
équivalent à B s’il existe deux matrices carrées inversibles Q' G C p et Q G 
C” telles que : 

A = QBQ'. 

Remarque 4.3.1 -Si A et B sont des matrices carrées équivalentes telles 
que 

A = Q~ 1 BQ 

les espaces Ai a et Ai b sont isomorphes, et pour tout polynôme. 

P(A) = Q~ l P(B)Q 


-Si P\ est le polynôme 


p i m = n (x - A'p 

i— 1 

où (Àj,..., A^,) sont des nombres complexes distincts, les espaces £p 0 et £p 1 
sont isomorphes. 

On veut définir maintenant / (A). Les fonctions envisagées posséderont 
la propriété (*) suivante : pour tout i tel que 1 < i < r, la fonction / 
possède des dérivées d’ordre m, — 1 en i. Par exemple : 

1. la fonction / est holomorphe dans un ouvert contenant les A* ; 

2. si les A i sont réels, la fonction / est dérivable rn„, — 1 fois en tout point 

A i , 

3. si les rrii sont tous égaux à 1, la fonction / est définie en A*. 
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Théorème 4.3.1 [59] S oit S une algèbre de fonctions possédant la propriété 
(*) et contenant les polynômes. Il existe un morphisme x d’algèbres de S 
dans Ai a, tel lue, pour tout polynôme P, on ait 

X (P) = P(A) 

Les valeurs propres de x (/) sont l es nombres f (Xi), associés aux mêmes 
sous-espaces propres que ceux de A (en prenant éventuellement la somme di¬ 
recte des sous-espaces propres associés à des valeurs propres ayant la même 
image par f). 


Preuve. Notons p l’application de £ dans C" définie par 

¥>(/)= (f W 0,)) 

c’est un morphisme d’algèbre et l’on pose 

x = 'ifoèp - 1 o Lp. 


On obtient ainsi un morphisme d’algèbre de S dans A4 a, et pour un poly¬ 
nôme P, on a bien 

X (P) = P(A). 

D’autre part, si 

X (f) = P(A), 

on a 

<p(f) =( P ( p ) , 

et donc 

f(Xi) = P(\i). 

La seconde propriété est vraie, puisqu’elle l’est pour les polynômes. 

On notera dans la suite 


X (f) = P{A). 


Remarque 4.3.2 -Ce qui précède permet de définir f (r) pour un endo¬ 
morphisme de matrice A. 

-A toute relation fonctionnelle correspond une relation matricielle ana¬ 
logue. Par exemple 

e tA e sA — e h+ s ) A 
ou 

sin 2 A + cos 2 A = I. 
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Quelques cas particuliers 

1. Si / est une fonction continue en A 

f(XI) = f(X)I. 

2. Si A à pour valeur propre unique A d’ordre rn. et si / est m — 1 fois 
dérivable en A, on a 


f<A) = f:CD(A-uy, 


p =o 


pi 


car le polynôme 


ml ,( p ) / y \ 

p ( x ) = E-A 1 ( x - x ) r ’ 


p=0 


pi 


est tel que, si 0 < j < m — 1; 


P U) ( A ) = fU) ( A ). 


3. Si les racines de P 0 sont toutes simples, on a 

f ( A \ = f ( \ ) ~~ M •" ~ ~ K+iI) ••• (^4 - A r I) 

J[ J tr ( A * - a o - ( A * - A -o ( A * - A *+i) - ( A * - A -) 

en partant du polynôme d’interpolation de Lagrange. 

4. Si A est inversible et si / (x) = -, on a / (A) = A -1 . 

x 

5. Si / = H est une fraction rationnelle n’ayant pas de pôles en A*, alors 

f(A) = P(A) .Q(A)- 1 . 

Théorème 4.3.2 [53] Soit f une fonction possédant la propriété (*) en A i 
et g une fonction possédant la propriété (*) en f (Ai). Alors g o f possède 
la propriété (*) en Xi et 


9° f(A) =g(f(A)). 
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Preuve. Il suffit de remarquer que si P est un polynôme tel que, pour i et 
j tels que 1 < i < r et 0 < j < — 1, on a 

f ü) (Xi) = P {j) (Xi) 

et que si Q est un polynôme tel que, pour i et j tels que 1 < i < r et 
0 < j < rrii — 1, on a 

9® (/ (Ai)) = Q (i) (/ (A,)) = <2® (P (Ai)) 

Alors 

(9°/)® (Ai) = (Q o Pf> (A,) 

et donc 

(g°f)(A) = (Q o P) (A) 

= Q(P(A))=g(f(A)) 


Corollaire 4.3.1 Si f admet une fonction réciproque f 1 avec les condi¬ 
tions de dérivabilité voulues, alors 

r i (f(A)) = A. 

Si les Xi sont positives et si 0 est au plus racine simple de / o, on peut 

, \ 2 

définir y/Â, |A| et l’on aura ( \/C4j = A et \TfË = |A| . 

Théorème 4.3.3 Soient A et B deux matrices semblables telles que A = 
Q~ 1 BQ, et f une fonction telle que, pour tout i compris entre 1 et r, f soit 
mi — 1 fois dérivable en Xi. Alors 

f(A) = Q~ 1 f(B)Q 
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Conclusion et perspectives 

Dans notre travail, on a donné le système des EDS des valeurs propres et 
des vecteurs propres de la solution d’une équation différentielle stochastique 
générale définie à l’aide d’un G—mouvement Brownien matriciel. Ce sys¬ 
tème a été difficile à obtenir à cause de la nature du G—mouvement Brow¬ 
nien et de la non-linéarité de la G—espérance. Ajoutée à cela, notre princi¬ 
pale difficulté réside dans le fait que les entrées d’une matrice G—Brownienne 
ne sont pas indépendantes en général. Pour contourner ces difficultés, on a 
supposé dans notre modèle que la variation quadratique (B) est une ma¬ 
trice diagonale. La G—formule d’intégration par parties matricielle a été la 
clé de ce travail. Un résultat intermédiaire sur le fait que les valeurs propres 
ne se rencontrent jamais a été également obtenu dans un cadre général 
Comme perspectives, nous envisageons de traiter le cas d’autres proces¬ 
sus liés à un G—mouvement Brownien et éventuellement considérer d’autres 
modèles. 
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